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Abstract — A modular method was suggested before to recover 
a band limited signal from tlie sample and hold and linearly 
interpolated (or, in general, an nth-order-hold) version of the 
regular samples. In this paper a novel approach for compensating 
the distortion of any interpolation based on modular method has 
been proposed. In this method the performance of the modular 
method is optimized by adding only some simply calculated 
coefficients. This approach causes drastic improvement in terms 
of SNRs with fewer modules compared to the classical modular 
method. Simulation results clearly confirm the improvement of 
the proposed method and also its superior robustness against 
I additive noise. 

I. Introduction 

Digital to analog converters are common in digital signal 
' processing and communication systems to reconstruct an ana- 
, log signal from its discrete time samples. Several methods with 
different names were introduced in the literature in 1970's and 
1980's f3l. S&H and LI were the dominant methods before 
that time; today. Polynomial interpolation and B-Spline are the 
] usual interpolation functions |i4l-||6|. 

These interpolators create some distortion at the Nyquist 
rate after low pass filtering, especially when S&H or LI are 
utilized. The advantage of these types of interpolators is their 
simplicity which makes them proper for practical use. To 
alleviate this problem, several methods such as inverse Sine 
I filtering, over-sampling, nonlinear and adaptive algorithms 
121-191, a modular method of the recovery of a signal from 
I its sampled-and-held version are described in [IJ for the 
■ uniform samples, 121 for the nonuniform samples, and succes- 
sive approximation using an iterative method fTOl-fT?! were 
introduced. The modular method is compared to the inverse 
Sine filtering in HI which shows that by using a few numbers 
of modules, the performance of the modular method excels the 
inverse filtering as far as noise is concerned. Over-sampling 
is not a practical solution due to its bandwidth requirements. 
The iterative method fTOl outperforms the modular method at 
the cost of more computation. 

We propose an Optimized modular method which enhances 
the performance of the classical modular method |1J. Our 
method is based on some optimum coefficients which are 
computed very simply by solving a least square problem. 
Indeed, these coefficients are calculated just one time for a 
specific interpolation system and are independent of signal 
to which the modular method is going to be applied. The 
coefficients themselves do not increase the complexity of the 



modular method and is very simple for practical usages. The 
simulation results show that the coefficients are well optimized 
and perform better then classical method. 

The rest of this paper is organized as follows: Section [III 
describes our general framework and introduces the terms and 
concepts used throughout the paper Section |lll] introduces our 
proposed method and is a straight forward manner to find 
the optimum coefficients for the modular method. Simulation 
results and comparison with the classical modular method for 
various interpolation systems will be presented in section HVl 
and finally, section |V] will conclude this paper. 

II. Preliminaries 

In this section we give a brief overview of the modular 
meffiod m that compensates the distortion of any interpolator 
such as Sample and Hold (S&H) and linear order hold by 
mixing the sum of cosine waves and then passing them through 
a lowpass filter. 

Suppose x(t) is sampled at the Nyquist rate (^) and assume 
s(t) is any interpolating function that fits the samples of 
x{t). According to these assumptions it can be formulated as 
follows: 

+ CX) 

sit) = hit) * ^ xinT)5it -riT) (1) 

«=-oo 

where hit) is the impulse response of the interpolation func- 
tion. The above equation can be written as shown below in 
frequency domain: 

Sif)^Hif)xY,Xif-i/T) (2) 

/=-oo 

According to the modular method, an improved reconstruc- 
tion of xit) can be drived from sit) by following process: 



m = ^(f) ^1 + Icosi—) + ■■■ + 2cosi^)j * n(f) (3) 

where n(f) is a lowpass filter with a bandwidth equal to 
the bandwidth of xit) iW). Eq. O) can be rewritten in the 
frequency domain as follows: 




Fig. 1. The reconstruction block diagram using Modular Method 
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= n(/) ^ H{f - j/T) ^ X(f - i/T - j/T) (4) 

J=-N /=-oo 

Since X{f) is band-limited (|4]l can be simplified as follows: 

+N 

X{f) = X(f) X n(/) H(f - j/T) (5) 

j=-N 

Therefore, it is obvious that as increases x will per- 
fectly converge to x if !!(/) Yj^^-n ~ -Z/^) becomes unity 
which always occurs for all interpolation functions since 
f^{Zp-c^H(f - j/T)} = 6(t), for instance this summation 
becomes summation of sine functions which is unity for 
Sample and Hold (S&H) interpolation function. However, 
practically it is not possible to apply infinite numbers of 
modules and only limited numbers of them are implemented 
by means of oscillators shown in Fig. [1] So the method will 
have distortion and to measure the distortion, let us define the 
mean-square error as: 



■ j/T) - 1 



df 



(6) 



and it is obvious that lim e - 0. 

Another issue with this approach is that the signals are 
considered to be analog. In practice, most of the signals that 
we are dealing with are discrete especially in a computer for 
interpolation of images and audio files. 

In the next section, we will drive all these relations again 
in the discrete domain and also minimize the mean square 
error using the optimization coefficients. The second part is 
the main part of this paper. 

III. Proposed Optimized Modular Method 

A. Modular method in Diecrete Domain 

Consider xln} is band-limited discrete signal which is down 
sampled and interpolated at the Nyquist rate {\/T) by h[n\ 



and the result is sin}. Therefore: 

s[n} = h[n} * x[kT}b[n - kT} (!) 

k=—oo 

and equivalently: 

+00 

S (k) = H(k) * Yj ^((^ - 'NfT))^ (8) 

/=— 00 

where S (k), H(k) and X(k) respectively are N-point DFTs of 
s[n}, h[n} and x[n}. Supposing that is divisible by T, the 
modular method can be formlated as follows: 



x[n} = LPF{s[«] (1 -t- 2cos(27m/T) + ■■■ + 2cos(2M7m/T)) } 

(9) 

where LPF is a FFT lowpass filter Also it can be shown that 
applying more than [T/2} modules not only does not enhance 
the performance but also may distort it. Consider T is even 
integer number and we have applied k modules more than T/2, 
then we have: 



i + i:J;'2co.(^) 



+k 



= 1 + Zli 2cos(^) + Z%i 2cos( 



1 + 2cos{^) + Zli (-l)"2cos(^) 



linn \ 



(10) 



It is abvious that the third term in the above final result will 
distort previous modules; this effect can be shown for odd Ts 
in the same way too. Hence, the maximum number of modules 
that is able to be applied is [T/2}. 

Furthermore, in the case that the maximum number of 
modules are applied, if we put the multiplicand of last cosine 
1 instead of 2, we will reach the impulse train. By means of 
Fourier series it would be proved as follows: 



JlnkniT 



(=-oo 
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T 



~6[n] = = 

X T 
T/2-1 

1 -H(-l)" + Y 2cos{2nin/T) 

1=1 

1 -H ^ 2cos{2nin/T) + cos{2n—n/T) 



(11) 



So, in this situation it will gather the original samples and if 
the filter is ideal, the output would be perfectly interpolatedQ 
The equation (|9]l will be rewritten in the frequency domain 
like below: 

+M 

X{k) = X{k) X n(^) H((k - jN/T))N (12) 

j=-M 

'Also true for analog D/A where 2,2, ...,2, 1 is not optimum but better 
than before. 



Simultaneously, the error of the interpolation process would 
be formulated in the same manner performed in the pervious 
section: 

+NIT ( +M 



-NIT 



e=Yu Z ^^^^ " jNIT))N - 1 



(13) 



Again minimizing this error is our main goal and is related 
directly to the performance of our system. In next subsection, 
our proposed Technique will be introduced to achieve this 
purpose. 

B. Optimization Coefficients 

Our goal is to minimize (fT3] l and as a result reduce 
distortion in the modular method in order to reach more 
precise interpolation for x. Our idea is that modules could be 
applied with some coefficients shown in Fig. |2] By choosing 
these coefficients appropriately, the performance of the method 
increases efficiently and we can achieve the same result with 
fewer modules. 

Now consider modules are multiplied by cf. 



2nn 2Mnn 
x[n] = LPF^I s[n] \ 1 + 2cicos{—) + ■■■ + 2cmcos( ^ ) 



and therefore, (flJl i becomes; 

+N/T f +M 



(14) 



-NIT \ j=-M 

For convenience the following parameters are defined: 



Hj(k) ^ H{{k-iNIT))N + H{{k + iNIT)% 
= F¥lN{h[n]x2cos{^)} 



' HiiO) H2(0) 
Hiil) H2(l) 

HiiN/T) H2{N/T) 



Hm(0) ^ 
Hm(1) 

HMiN/T), 



(16) 



(17) 
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Fig. 2. The reconstruction block diagram using Discrete Modular Method 
and Optimazation coefficients 
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Now from these new definitions the error function can be 
simply rewritten as follows: 

e = (norm(§C; - 93))^ (19) 
Hence, to minimize this error, we should solve: 

§(£ = 93 (20) 

Since the number of equations (N/T) are much more than 
the unknowns (M), we could find the optimum answer by 
considering the mean square error method and the pseudo- 
inverse: 



(21) 



Therefore, the problem can be solved and we have found 
some coefficients which minimize the error for the finite 
number of modules without any exception on interpolation 
function. Furthermore, we will show in next section that 
these coefficients cuase dramatic results in comparison to 
the classical modular method which is a special case of our 
method by assigning the coefficients one. The main key in our 
method is that the coefficients are calculated very easy and 
fast and by only accessing impulse response of interpolation 
function (h[n]). Moreover the coefficients are calculated just 
one time for an interpolation function and are stored in a 
lookup table and does not need to find them again every time 
we need them. 

IV. Simulation Results and Discussion 

We utilized MATLAB® simulation environment to evaluate 
and compare the performance of methods. To have fair com- 
parison, initial band limited signals are produced randomly, 
and the performance of each method is averaged over 100 
signals. The initial signal is FFT lowpass filtered version of 
white Gaussian noise signals. To show the significance of this 
method, the sampling rate is performed at the Nyquist rate. 
The performance criterion for our simulations is the Signal to 
Noise Ratio (SNR) in dB. To avoid transient errors at the 
end points, SNR is calculated for interior points and 10% 
of the end points are ignored. As illustrated in Fig. |3] the 
SNR increases monotonically in dB for classical method as 
the number of modules increases, while the optimum method 
increases exponentially as the number of modules increases. 
This means more than 250dB for simple S&H interpolation 
and this is quite impressive in real engineering applications. 

Fig.Hlshows similar results for the Linear Interpolation (LI). 
The difference between the classical method and the optimum 
at the first few numbers of modules is not very significant. 
However, as the number of modules increases, the difference 
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Fig. 3. SNR vs. the number of modules for Classical Modular method and Fig. 4. SNR vs. the number of modules for Classical Modular method and 
our proposed Optimized Modular method for S&H interpolation our proposed Optimized Modular method for Linear interpolation 



becomes apparent. This shows that our method can find the 
optimum coefficients for any interpolation independently. 

To study the eff'ect of noise, we added a white Gaussian 
noise to the band limited signal. This is the model of the 
electronic devices that generate thermal noise. Fig. |5] shows 
that for low SNR, the performance of our method is not so 
significant over classical method, however as much as the 
power of noise decrease our method give much greater SNRs. 

These simulations show that using the modular method 
by means of optimized coefficients enhance its performance 
dramatically boffi in noisy and noiseless environments inde- 
pendent of ffie type of interpolation. 

V. Conclusion 

A novel Optimized Modular method is proposed for com- 
pensating error of any interpolation system. We add the 
optimized coefficients calculated in a very simple manner into 
the Classical Modular method in order to maximize its perfor- 
mance as much as it could be. The simulations show very high 
improvements versus classical methods. This proposed method 
is also more favorable in terms of computational complexity 
with respects to other error compensating methods, since not 
only the modular method has very simple algorithm but also 
we could reach better SNR values with just 2 modules rather 
than 5 modules in the classical method. 
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